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Abstract. 

For uniaxial easy axis films, properties of magnetic domains are usually 
described within the Kittel model, which assumes that domain walls are much 
thinner than the domains. In this work we present a simple model that includes 
a proper description of the magnetostatic energy of domains and domain walls 
and also takes into account the interaction between both surfaces of the film. Our 
model describes the behavior of domain and wall widths as a function of film 
thickness, and is especially well suited for the strong stripe phase. We prove the 
existence of a critical value of magneto-crystalline anisotropy above which stripe 
domains exist for any film thickness and justify our model by comparison with 
exact results. The model is in good agreement with experimental data for hep 
cobalt. 
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1. Introduction 

Thin magnetic films with stripe domains and perpendicular magneto-crystalline 
anisotropy present a high scientific interest. They are model systems to understand the 
domain structure of ferromagnetic, as well as ferroelectric, thin films. Such materials 
are used for the fabrication of memories and spin injection devices. Considering a thin 
film with its bulk easy axis perpendicular to the surface, the magnetization direction 
in the domains can depend on the film thickness. Thick films behave as bulk material 
but for thin films, the competition between exchange, magnetostatic and magneto- 
crystalline energy can tip up the magnetization into the plane. A very well known 
theory of magnetic domain structure has been developed by Kittel jT]. In this model, 
the film magnetization is described by a square profile. Domain wall energy is allowed 
for the energy balance, but not for the magnetization profile, which determines the 
magnetostatic energy. Such an approximation is only valid when the domain width is 
much larger than the width of domain walls. Therefore, this model is bound to fail 
for very thin films, where the domain width vanishes. Furthermore, the magnetostatic 
interaction between the top and bottom surfaces is disregarded, which leads to wrong 
results unless the domain width is much less than the film thickness. The latter does 
not hold for very thin strongly anisotropic films. 

Within the Kittel model, the width of stripe domains is proportional to the square 
root of the film thickness. The model also predicts the existence of a critical film 
thickness below which the magnetization direction flips from out-of-plane to in-plane 
direction. Later on, it was shown that for materials with strong magneto-crystalline 
anisotropy K (in comparison to the square of the saturation magnetization Mq) 
the magnetostatic interaction energy between top and bottom interface cannot be 
neglected. The ratio 2K / {^qMq) (in SI units) is also known as the quality factor Q 
and the magnetization direction remains perpendicular to the surface for all values of 
film thickness for Q > 1. Numerous theoretical and experimental works have been 
devoted to such materials with strong perpendicular magnetic anisotropy [21 [3j |4] . It 
has been shown that in such materials the domain width has a minimum as a function 
of film thickness, increasing towards smaller or larger thicknesses. Materials with 
strong perpendicular anisotropy are also of technological importance in the field of 
magnetic multilayers (see W. and references therein). 

Just above the critical film thickness, the magnetization still lies predominantly 
in the plane of the film, deviations from the homogeneous in-plane orientation being 
very small. This permits an exact mathematical solution of the problem of stripe 
domain nucleation at the critical point [B]. Near the critical point the out-of-plane 
component of the magnetization remains much smaller than the saturation value Mq 
and there is a strong variation of the magnetization direction across the film thickness. 
The corresponding phase is usually called the weak stripe phase. As the film thickness 
increases, the magnetization direction in the domains tends more and more towards the 
bulk easy axis, which is perpendicular to the film. Gradually, the weak stripe phase 
evolves towards the strong stripe phase, where the magnetization magnetization is 
predominantly parallel (or antiparallel) to the bulk easy axis. 

For a detailed analysis of the weak stripe phase a micromagnetic numerical 
analysis is indispensable (see [3 El El [10] and many other works). However, simple 
analytical models are helpful for getting a quick, albeit crude understanding of the 
stripe phases. These models restrict the magnetization direction to the vertical 
plane parallel to the stripes (i.e. they consider domain walls of Bloch type) but 
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improve the Kittel model. One can distinguish between models with a linear domain 
wall profile [TT] or others which use the Jacobi sine function to parametrize the 
magnetization profile [TH [13] . Also a sinewave magnetization profile was proposed 
by Saito [14, to treat the weak stripe phase in an approximate manner. In the 
following we present a thorough analytical and numerical analysis of a simple but 
quite complete model for the strong stripe phase. We adopt a sinewave magnetization 
profile in the wall and we take into account the magnetostatic interaction between 
the top and bottom surfaces of the film. Our model applies to hexagonal Co, for 
which a critical thickness between 25 nm (TSl US] and 40 nm [T7] has been reported, 
as well as to materials with high magneto-crystalline anisotropy, such as FePd(OOl) 
or garnet films [U [T7J [18] . Another possible application concerns ferromagnetic thin 
films of MnsGes that were recently synthesized [TU[2n]- Our model unifies previous 
results for materials with strong anisotropy {Q > 1) [U[3J|1] and those with medium 
anisotropy (5 ^ Q < 1) [H], containing them as special cases. The explicit treatment 
of the domain wall is an important improvement with respect to the Kittel model 
since it allows the determination of the in-plane remanent magnetization, which is 
important for the analysis of experimental data. We give numerical results for the 
dependence of domain and wall widths on the film thickness and calculate the critical 
thickness at which the magnetization switches from out-of-plane to in-plane direction 
if Q < 1. We show that the critical thickness tends to zero as Q — ?> 1. This behavior 
is identical in our model and in the exact solution [6] . In the case of strong anisotropy 
(Q > 1) stripe magnetic domains exist at any film thickness. 

Our theory considers a non-magnetistrictive film: the magneto-elastic energy is 
equal to zero. Moreover, the surface anisotropy is neglected. It is not suited for 
materials where the surface anisotropy plays an important role. We concentrate on 
the strong stripe phase since the details of the weak stripe phase are hard to capture 
in a simple model. In materials with high anisotropy the weak stripe phase does 
not appear at all and it will be shown below that even in materials with medium 
anisotropy < Q < 1) it is of rather limited importance. In the latter case the weak 
stripe phase occurs only in a narrow interval of thicknesses around the critical point. 
This region of the weak stripe phase is not correctly described in our model. On the 
other hand, even in the case of weak anisotropy {Q ^ 1), our model provides a correct 
qualitative description of the strong stripe phase, although it neglects such important 
phenomena as closure domains. 

Our paper is organized as follows. After a short summary of known results 
(Kittel's theory, sawtooth magnetization model, exact results for stripe phase 
nucleation) we present our model in Section 3. In Section 4 we show its analytical 
and numerical solution. The usefulness of our theory is demonstrated in Section 5 
by way of theoretical analysis of published experimental data for the well-known hep 
cobalt system. In Section 6 we compare our model with other approaches and expose 
its strengths and limitations. 

2. Known results 

2.1. Kittel's model 

This model was developed for ferromagnetic films with uniaxial anisotropy 
perpendicular to the film. It determines the width of magnetic stripe domains d as 
a function of the film thickness h, and the critical thickness where the magnetization 
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Figure 1: Kittel model for ferromagnetic thin films with uniaxial anisotropy. 



unit volume) is given by: 

£;totai = 0. 136^0 Afo'^ + ^, (1) 

where Mq denotes the saturation magnetization and a the surface wall energy. The 
latter is obtained in a separate variational calculation for an isolated Bloch wall [H] , 
a = 4\/ AK, where A is the exchange stiffness constant. The same calculation yields 
the wall width, S ~ tt-^/ A/ K. It is convenient to express all spatial dimensions in the 
units of y^A/K: 

h-^^, d-^^, ^ = ^=- (2) 
y/A/K' ^A/K' ^A/K 

Thus, Kittel's reduced wall width is a constant, 

5 = TT. (3) 

The only essential material parameter in the model is the dimensionless quality ratio, 
Q = 2K/{^oMq). The first term of Eq. (jlj refers to the magnetostatic energy of 
a rectangular domain profile, neglecting the influences of finite wall width and film 
thickness. Therefore, strictly speaking, the Kittel model is only applicable when the 
following strong inequalities hold: 

S <^d<^h, (4) 

or alternatively, S <ti d ^ h. The right-hand condition ensures the negligibility of the 
magnetostatic interaction between the top and the bottom surfaces of the film. By 
minimizing the energy, one obtains the equilibrium domain width: 

d = SMy/ol, (5) 
and the minimum energy density: 



= 2 f''Y'o^ . (6) 

As one can see, the minimum energy density E^™j tends to infinity as /i — > 0. 
Therefore, the multi-domain structure cannot be stable in that limit within the Kittel 
model. Indeed, for a uniform in-plane magnetization (perpendicular to the bulk 
easy axis), the magnetostatic and the exchange energies are equal to zero, and the 
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anisotropy energy is maximal and equal to K. The critical point is obtained by 
equating Eq. ^ to K . Therefore, we can express the reduced critical thickness as 
follows: 



Within the Kittel model, the domain width is identical to the half-period d (see Fig[T]). 
Its value at the critical point is obtained by setting Eq. ([7]) into dU, which results in: 



These equations permit to obtain a qualitative evaluation of the critical quantities. 
But they fail in the neighborhood of the critical thickness, where domain and wall 
widths become close, cf. Eqs. ([3]) and ([8|). Furthermore, by Eqs. ([7]) and dU in the 
critical region the ratio h/d equals 0.544/Q — 0.272/ioMg /if , i.e. also the second 
precondition of the Kittel model (jl} is violated in the case of large anisotropy or small 
magnetization, Q ^ 1. As a result, Kittel's model makes a wrong prediction for large 
Q: according to Eq. ([7]), ^" decreases asymptotically, remaining nonzero at any finite 
Q. In actual fact, the critical thickness should vanish at Q = 1, as shown by rigorous 
calculations [6]. 

For that reason we decided to develop a model that would not break down in 
the strong-anisotropy case, by allowing for the electrostatic interaction between the 
top and bottom surfaces of the film and taking into account the wall width explicitly. 
Within that improved model we have obtained a good description of domain and wall 
width variation as a function of the film thickness (see Section 3). 

2.2. Sawtooth magnetization model 

The main goal of Kittel's model is the description of the strong stripe phase. Yet there 
is no direct transition between strong stripes and planar magnetization. Rather, it 
proceeds via an intermediate weak stripe configuration. The latter can be described 
by sawtooth magnetization model [TJ. In this model the canting angle of the 
magnetization out of the homogeneous in-plane direction, 6(x), is assumed to be a 
saw-tooth curve with a maximum value of 0o and a minimum value of — Gq, 6o < 7r/2. 
Between the extrema the variation is linear and the half period is denoted by d. The 
sawtooth magnetization model is an approximate one, since it neglects the variation of 
the magnetization direction across the film thickness. The magnetostatic interaction 
between the top and bottom surfaces is not taken into account either. From the 
equation for the total energy below one can obtain the expression for the critical 
thickness where the magnetization starts to fiip out-of-plane by a small deflection 
angle. The first part of the total energy corresponds to the magnetostatic energy; 
Co is the leading coefficient of the Fourier series obtained from periodic profile of the 
deflection angle. The maximum angle 0o is a variational parameter, < < 7r/2. 




4.35Q 



-I 



(7) 




(8) 




(9) 



The critical thickness is expressed as : 



(10) 
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The expression above has the same form as Kittel's critical thickness ft.", Eq. (O, but 
with a smaher pre-factor. This opens up a possibihty to interpret /i" and /i"' as the 
upper and lower bounds of an interval where the weak stripe phase exists. The latter 
is a transitional state between the homogeneously magnetized in-plane configuration, 
stable below ft"', and Kittel's strong stripe domain structure, taking place above ft". 
One observes the 'wrong' behavior of h'^'{Q) at Q > 1, similar to that of h'^{Q). Eq. 
(|10|) also misbehaves at Q <C 1. Rigorous calculations show (see Section 2.3) that the 
true lower bound of the weak stripe phase does not diverge as Q 0. 

2.3. Exact description of stripe domain nucleation 

The exact solution of the nucleation problem was given in 1961 by Muller [6]. A 
modern presentation, summarized here, can be found in the book of Hubert and 
Schafer [7^. The theory was formulated for a thin magnetic film (thickness ft, 
magnetization Mq) with uniaxial anisotropy perpendicular to the film (anisotropy 
constant A'). The energy expression contains the magnetostatic energy, the exchange 
(stiffness A) and anisotropy contributions; the theory imposes no restrictions on the 
magnetization direction. Below the critical thickness he the magnetostatic energy 
forces the magnetization into the plane of the film. At h ^ he there is an instability and 
the weak stripe phase emerges. In the neighbourhood of the critical point, deviations 
from a homogeneous in-plane orientation of magnetization are small, which allows to 
linearize the system of micromagnetic equations and find the exact analytical form of 
the instability mode. The reduced critical thickness he is a universal function of the 
quality ratio Q, see Fig. |4] (solid line). (The quantity plotted in the original drawing. 
Fig. 3.109a of Ref. [7], is a factor 2tt smaller.) At the same time, the critical half- 
period dc can be found, as well as the distribution of magnetization directions in the 
critical mode. In the weak-anisotropy limit, Q — ?> 0, both he and dc tend to the same 
finite value, 

he de ^ 2tt. (11) 

As Q increases, the critical thickness he decreases and vanishes at Q = 1, whereas de 
diverges at that point. The exact results will be used later on to judge the validity of 
the sinewave wall model. 

3. Sinewave wall model (SWM) 

We propose a model stripe domain structure as shown in Figs. [2]and[3l Like in Kittel's 
model, the half-period is denoted by d, but it now includes a domain wall of width 
6 (see Fig. [5]). The inner domains, with constant magnetization equal to ±Mo, have 
a width oi d — S. The walls separating different domains are assumed to be of Bloch 
type, with linear variation of the angle Q{x) between the magnetization direction 
and the y-axis. This results in a magnetization profile of sinewave form, which is a 
good approximation of the profile obtained by the variational method. The schematic 
representation in Fig. [5] shows that we consider parallel stripe domains infinite in the 
y direction. In the x direction, the periodic magnetization profile is as shown in Fig. [3] 
The z dimension is restricted to the film thickness ft. 

To calculate the magnetostatic energy, we use an analogy with the electrostatic 
field calculation for alternating, positively and negatively, charged stripes (see Landau- 
Lifschitz [22]). We consider the realistic magnetization profile of Fig. [3] (as opposed 
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to a simple rectangular meander) as well as the magnetostatic interaction between the 
top and bottom surfaces of the film. 




Figure 2: Sketch of periodic stripe domains in a magnetic thin film with domain walls 
of finite thickness. 
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Figure 3: (color online) One full period of the domain structure of Fig. [2] M^^ix) 
(solid line) denotes the z component of the film magnetization, and Q(x) is the angle 
between the magnetization direction and the film surface (dashed line). The domain 
wall width is denoted by 8. 



The magnetostatic energy density (per volume) can be written as: 
fe=i 



1 — exp — Trfc 



with Ck 



2Mn 



ki: [l-k^iS/d) 



21 ^ 2d 



d 
knS 



(12) 
(13) 
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where k is an odd number, Ck is the Fourier coefficient of Mz(x) (see Fig. [3]). The 
domain structure is determined by the interplay between magnetostatic, exchange and 
anisotropy energies. The latter two contributions are as follows: 



1 f'" , /de^ ' 



- dS^ 



1 r'^'' 

i^aniso =^J iCC0s2[e]d; 



(14) 



(15) 



The half-period d and the wall width S are then calculated by minimizing the total 
energy, 

^-total = -E-mag + ^-cxch + -E-aniso ■ (16) 

Since it is not possible to solve analytically the two variational equations 9d£'totai = 0, 
dsEtotsd = in the general case, we use a Nelder-Mead numerical method. The results 
are presented in the next Section. 



4. Analytical and numerical results 

4.1. Large film thickness 

Before presenting the numerical solution in the general case let us discuss some specific 
situations. For sufficiently thick films the strong inequalities (HJ are fulfilled and the 
expression for the Fourier coefficients, Eq. (|13p. becomes: 



Then Kittel's form is recovered for the magnetostatic energy density: 

Em.g = ^°3^° ^ ^ 0.136moM2-. (18) 
fc=i 

The wall energy associated with the above magnetostatic expression is a minimum for 
the following reduced wall width: 

~5oo = 7rV2. (19) 

The corresponding surface wall energy, a — ttv^a/AK, is slightly (11%) larger 
than the variational result, 4\/ AK, because of the imposed sinewave wall profile of 
Mz{x). In the thick-film limit the dependence of the domain width on film thickness 
can be expressed analytically, 

(J=4.05y^, (20) 
and is very similar to the result of the Kittel model, Eq. ([5]) . 
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4-2. Critical thickness and critical anisotropy 

As the film thickness is reduced, the domain width d decreases and the waU width 
6 increases. Just above tlie critical thickness, the magnetization profile is purely 
sinewave (S — ?> d), forming an unstable spiral magnetic configuration. In this limit 
Eq. ([T2l) can be simplified: the sum disappears since only the leading term of the 
Fourier series, with k — I, survives. The magnetostatic energy becomes: 

(21) 



1 ( ^ 

1 — exp ( ^'"'^ 



and anisotropy plus exchange terms (|14I15I) . can be written as: 

^ K n'^A ^ ^ 

-Ewall - Y + (22) 

Minimizing the total energy, i?niag + £'waii, with respect to d and equating the result 
to K yields the critical thickness h'^: for this film thickness the total energies of the 
stripe domain structure and of a mono-domain state with in-plane magnetization are 
equal. The corresponding value of the half-period is dj,. Both /ij, and d'^ depend on 
the quality ratio, Q — 2K/{^oMq). This dependence can be presented in parametric 
form by introducing an auxiliary quantity, t = irh'^/d'^. One then finds 

(23) 



3 


'1 






2 


t 







~ r- 1-(1 + V3)e-* 
^-^^^V l-(l + t)e-* ' 

and 

4 = TTi-i/i^. ^ (25) 

The parameter t runs from zero to infinity. The resulting H'^-ys-Q dependence is 
displayed in Fig. |4] (dashed curve). It can be regarded as a borderline between the 
strong- and the weak-stripe phases. The solid line in Fig. |4] presents the exact result 
of MuUer's theory [SI 13: hc{Q). This line separates the region of weak stripes from 
the area of homogeneous in-plane magnetization, as observed in very thin films with 
Q < 1. 

For Q small, Q — > 0, t — > oo, the critical thickness h'^ diverges, 

K ^ 3^^Q-^ K 3.67Q-\ (26) 

whereas 

4 irVe. (27) 
The opposite limiting case is i — > 0, Q — > 1. In that limit h'^ tends to zero. 



h'^ = 4^2(1 -Q) , (28) 
while d'^ diverges. 



4 = -Vt^- ^''^ 

For strong anisotropy, Q > I or K > ^^^Mq, there is no physically meaningful h'^ or d'^ 
and, without a magnetic field, the mono-domain structure with in-plane magnetization 
is unstable for any h. We find it remarkable that our model reproduces the exact value 
for the critical quality ratio, Q = 1 [Hll^. The Kittel model is limited to systems with 
small quality ratios. However, in many strongly anisotropic materials Q is large and 
no critical thickness is observed [H [181 [23] . 
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Figure 4: Phase diagram of thin magnetic films in zero magnetic field. The solid 
line is the critical thickness he obtained in the exact approach [H |7] , it separates the 
homogeneous state without domains from the weak stripe phase. The dashed line is 
the reduced critical thickness h'^ of SWM, computed using Eqs. (^5]) and ([M)) . This 
is a cross-over line between the weak and the strong stripe phases. 



4-3. Numerical results 

Away from the critical region, the total energy (jl6p has to be minimized with respect 
to d and S. In general, for arbitrary Mq, K, A, and /i, the result cannot be expressed 
analytically and we have to resort to a numerical procedure. The dependence on A and 
K can still be taken into account by changing over to the dimensionless variables ([2|). 
Then we are left with two quantities 5 and d—5 versus h with only one parameter Q. In 
Fig. [5] we prefer to plot the difference d~~ 5 rather than d since d — 5 can be regarded 
as an order parameter. The numerical curves confirm all analytical expressions, in 
particular, those describing the asymptotic behavior at ^ 1, Eqs. (|19p and ([20)) . 
Thus, 5{h) can be seen to tend to a universal limit. 

One observes in Fig. [5] two distinct regimes, for Q < 1 and for Q > 1. If Q > 1, 
the stripe domain structure is stable at any h. The half-period of the structure, d{h) 
or has a minimum at a certain finite thickness, increasing towards smaller and 

larger h. A prominent feature of the curves with Q < 1 is the presence of a critical 
thickness h'^ where the width of the inner domains, d — 5, vanishes. At that point the 
sample contains nothing but domain walls, the magnetic structure being a spiral of 
period 2dc = 25c- Obviously, the wall width cannot be neglected, especially at Q ss 1. 
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Figure 5: Inner domain width d — 6 and wall width 5 as a function of film thickness 
h. All figures are presented in reduced coordinates ([2]). Bold numbers correspond to 
the value of Q. 



4.4- In-plane remanent magnetization 

In a magnetic field directed in the plane of the film the stripe domains arrange 
themselves parallel to the external field (see Refs. [HI US]). Each domain wall 
has a magnetic moment component parallel to the in-plane field direction, all of 
them being directed in the same sense for a small but finite external field. These 
contributions sum up to a finite in-plane remanent magnetization with a corresponding 
hysteresis, as observed in thin hep Co or MnsGes films [20l[l6]. This remanent in- 
plane magnetization (when the external field tends to zero) is easy to calculate in our 
model. The sinewave magnetization profile in the domain wall leads to the following 
expression: 



At the critical thickness, the remanent magnetization is equal to 2/7r ~ 63.7 % of 
the saturation magnetization. Just below this point we would expect the weak stripe 
phase [14] in a narrow interval of thicknesses before the transition into the planar 
mono-domain state. In the weak stripe phase the maximum z-component of the 
magnetization is less than Mq, like in sawtooth magnetization model [14) . However, 
the weak stripe phase is beyond the scope of our model. 




(30) 
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5. Hexagonal Co 

Cobalt hexagonal thin films have been intensively studied. Below a critical thickness, 
such films exhibit planar magnetization. Recent experimental publications [ITl [161 El] 
report critical thicknesses between 25 and 50 nanometers. Since the study of 
Brandenburg et al. [17] is the most complete one, SWM has been tested using those 
data, especially the half-period as a function of film thickness, shown in Fig. [6] 




Figure 6: Half-period d (full line) and wall width 6 (dashed line) as functions of 
thickness for Co films. Red points are experimental d values from Ref. [ITj. Red and 
black curves are fits to the Kittel model and SWM, respectively. 



SWM has three parameters to characterize different materials: the magneto- 
crystalline anisotropy constant K, the saturation magnetization Mq and the exchange 
stiffness constant A. The saturation magnetization was fixed to the experimental 
value of Mq = 1.43 MA/m, while the anisotropy as well as the exchange constant 
were chosen as fitting parameters. The best fit was obtained with K = 820 kJ/m'^, 
which is rather close to the bulk anisotropy constant (see Tab.[T]), and A = 45 pJ/m, 
which exceeds the value deduced from inelastic neutron scattering experiments |26j . 
However, it should be taken into account that A is not known accurately, since relating 
neutron scattering data with our continuum model involves several approximations. 

For comparison, we also fitted the experimental data to the Kittel model (see 
Fig. [SJ, leading to the parameter values presented in Table [Tj The saturation 
magnetization was fixed to the same value, but the anisotropy and exchange constants 
deviate more strongly from the experiment than those obtained in the SWM fit. This 
suggests that the Kittel model is less realistic than SWM. We also observe important 
differences between both models near the critical thickness. The curvature of d{h) 
near he is different in the two models. In SWM, for h between 40 and 80 nm the wall 
is wider than the inner domain, whereas in the Kittel model the wall width remains 



Theory of magnetic domains in uniaxial thin films 



13 



cobalt 


Kittel 


SWM 


measurements 


K [kJ/m3] 


1200 


820 


820", 920'' 


A [pJ/m] 


76 


45 


28" 


Mo [MA/m] 


1.43* 


1.43* 


1.43"* 


he [nm] 


36 


36 


36 ±3 


dc [nm] 


63 


67 


63 ±5 



Table 1: Summary of important physical constants obtained for cobalt. The asterisks 
indicate fixed parameters. Experimental data are extracted from: " Ref. [27], Ref. 
[15], " Ref. [26], Ref. [Ill [25] 



constant and always inferior to the inner domain width. 

The quality factor for Co is Q = 0.64, which is less than one. We expect that for 
materials with stronger anisotropy, i.e. with larger Q, the difference between SWM 
and the Kittel model will be even more significant. 

6. Comparison 

Our model is constructed in such a way that it applies at any Q, including the strong- 
anisotropy case, Q > 1, where there is no critical thickness. 

The most interesting region is the one of medium anisotropy, ^ ^ Q < 1- To 
illustrate that region, we chose hexagonal cobalt as an example. With the parameters 
of Table 1 (SWM) we evaluated the critical thickness in SWM, h'^ = 36 nm. This 
agrees with the experiment of Brandenburg et al. [T7] , who find a critical thickness of 
about 40 nm. However, in Muller's exact theory [B] the critical thickness comes out 
much smaller, = 23 nm. One has to take into account, though, that stripe domains 
were observed by other groups in Co films as thin as 25 nm and that weak stripes are 
certainly hard to see. Our interpretation is that the interval between Muller's he = 23 
nm and the SWM h'^ = 36 nm is a region of the weak stripe phase. 

In the case of cobalt, such an interpretation is supported by a calculation of the 
weak stripe phase within the sawtooth magnetization model. The two variational 
parameters are the maximum canting angle Oq and the half-period d. For better 
precision, we included the complete Fourier series as well as the magnetostatic 
interaction between the top and bottom surfaces. The numerical results are presented 
in Fig. [7J The weak stripe phase sets in at /i = 28 nm and has a lower energy than the 
strong stripe phase up to 44 nm. Between 28 nm and 44 nm the maximum canting 
angle increases from zero to 76 degrees. So, in the case of hexagonal cobalt, wc would 
expect the existence of the weak stripe phase between 23 and about 44 nm. At larger 
thicknesses the strong stripe phase prevails and SWM is more appropriate than the 
sawtooth magnetization model. 

Let us now briefly discuss the region of very small anisotropy, Q <C 1 . It allows 
mathematical simplifications, as already discussed in Section 4. There, it was shown 
that the range of thicknesses where the weak stripe phase could be expected becomes 
wider according as Q decreases, i.e. as the magnetic anisotropy weakens. For very 
weak anisotropy we also expect other complications, such as the presence of closure 
domains. We believe our model can still be applied if the layer of the closure domains 
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Figure 7: Comparison of the total energy of the weak stripe phase (sawtooth 
magnetization model, solid line) with that one of the strong stripe phase (SWM, 
dashed line). 



is much thinner than the film as a whole. 



7. Conclusion 

We present a model (SWM) for the strong stripe phase in magnetic thin films with 
perpendicular anisotropy, which improves the well-known Kittel model in two respects. 
Firstly, our model includes the magnetostatic interaction between the top and bottom 
surfaces of the film. This is important for materials with strong magnetocrystalline 
anisotropy. Secondly, the domain walls are treated explicitly, which improves the 
numerical accuracy and allows the calculation of the in-plane remanent magnetization. 
Our model is simple, but general enough to permit a thorough analysis of the strong 
stripe phase. It is especially suited to interpret experimental data. We demonstrate 
the existence of a critical anisotropy above which stripe domains exist at any film 
thickness. SWM reproduces the exact threshold value, Q — 1. We derive a number of 
analytical results facilitatng the estimation of important parameters. The numerical 
results obtained for different quality factors Q show the evolution of domain and wall 
widths as functions of film thickness. 

With the new model we are able to correctly describe the behavior of magnetic 
domains in cobalt thin films. It allows to explain the experimental observations of 
Brandenburg et al |17| with more realistic fitting parameters than using the Kittel 
model. By comparing the exact critical thickness for stripe nucleation he with our 
critical thickness h'^, which corresponds to the on-set of strong stripes, we are able to 
estimate the range of thicknesses where weak stripes are expected. We show that for 
Co this interval is rather narrow, 10 to 15 nm, which validates our model. We should 
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remark that the currently available experimental data for Co do not suffice for locating 
the thickness range of the weak stripe phase accurately. We find it important that 
our model takes the wall width into account in an adequate way. In the neighborhood 
of the critical thickness, domain and wall width are of comparable size. Thanks to 
its universality, our model can be applied to other types of ferromagnetic films, e.g., 
FePd or MusGea. 
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